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13.  AtSTIUCT  {MiMiitMniOOiMPdtf  This  document  represents  a  final  technical  report  for 

AFOSR  Grant  No.  89-0185.  The  research  program  supported  under  this  grant  involved 
a  fundamental  study  of  the  mechanical  properties  of  microelectronic  thin  film 
materials.  The  focus  of  the  work  was  on  the  microscopic  processes  that  lead  to 
stresses  in  microelectronic  thin  films  and  control  the  mechanical  properties  of 
these  materials.  Our  work  ranged  form  studies  of  interconnect  metals,  passivation 
glasses  and  heteroepitaxial  thin  film  semiconductors.  In  previous  reports  we 
highlighted  work  done  of  the  mechanical  properties  of  passivation  glasses  and 
interconnect  metals  and  the  processes  of  failure  of  interconnect  metals.  In  this 
final  report  we  highlighted  our  work  on  misfit  dislocation  formation  in  hetero- 
pitaxial  thin  films.  The  body  of  the  report  contains  two  reports  of  our  work  on 
misfit  dislocations.  Our  overall  understanding  of  misfit  dislocations  in  Si-Ge 
layers  is  described  in  the  first  report.  That  report  also  descirbes  our  work  on  the 
effect  of  capping  layers  on  the  formation  of  misfit  dislocations.  The  second 
report  descirbes  our  in  situ  studies  of  the  kinetics  of  formation  of  misfit 
dislocations  in  Si-Ge  films.  By  measureing  the  change  of  substrates  curvature  as  a 
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function  of  time  during  anneAling  we  were  able  to  determine  the  way  in  which 
mobile  threading  dislocation  density  changes  during  the  course  of  annealing. 
This  work  has  led  to  a  much  more  complete  understanding  of  the  kinetics 
of  strain  relaxation  in  heteropitaxial  thin  films. 

The  presentations  and  publications  that  have  resulted  from  this  work  are  listed 
at  the  end  of  this  report. 
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I.  SUMMARY 


This  document  represents  a  final  technical  report  for  AFOSR  Gram  No.  89-0 1 1 85.  The  research 
program  supported  under  this  grant  involved  a  fundamental  study  of  the  mechanical  properties  of 
microelectronic  thin  film  materials.  The  focus  of  the  work  was  on  the  microscopic  processes  that 
lead  to  stresses  in  microelectronic  thin  films  and  control  the  mechanical  properties  of  these 
materials.  Our  work  ranged  from  studies  of  interconnect  metals,  passivation  glasses  and 
heteroepitaxial  thin  film  semiconductors.  In  previous  reports  we  highlighted  work  done  of  the 
mechanical  properties  of  passivation  glasses  and  interconnect  metals  and  the  processes  of  failure  of 
interconnect  metals.  In  this  final  report  we  highlight  our  work  on  misfit  dislocancai  formation  in 
heteroepitaxial  thin  films 

The  body  of  this  report  contains  two  reports  of  our  work  on  nrisfit  dislocations.  Our  overall 
understanding  of  misfit  dislocations  in  Si-Ge  layers  is  described  in  the  first  report  That  repon  also 
describes  our  wcn-k  on  the  effect  of  capping  layers  on  the  formation  of  misfit  dislocations.  The 
second  repon  describes  our  in  situ  studies  of  the  kinetics  of  fmmaticn  of  misfit  dislocations  in  Si- 
Ge  films.  By  measuring  the  change  of  substrate  curvature  as  a  function  of  time  during  annealing 
we  were  able  to  determine  the  way  in  which  mobile  threading  dislocation  density  changes  during 
the  course  of  annealing.  This  work  has  led  to  a  much  more  cocqplete  understanding  of  the  kinetics 
of  strain  relaxation  in  heteroepitaxial  thin  films. 

The  presentations  and  publications  that  have  resulted  from  this  work  are  listed  at  the  end  of  this 
report. 


II.  RESEARCH  REPORT 

A.  Mechanisms  and  Kinetics  of  Misfit  Dislocation  Formation  in 
Heteroepitaxial  Thin  Films 
W.D.  Nix.  D.B.  Noble  and  J.F.  Turlo 

1.  Introduction 

A  great  amount  of  attention  has  been  focussed  recendy  on  the  processes  by  which  misfit 
dislocations  form  in  heteroepitaxial  structures;  see,  for  example,  [1-2].  A  number  of 
semiconductor  device  applications  tely  on  the  elastic  strain  in  heteroepitaxial  films  to  control  the 
bandgap  in  these  materials  [3-4].  For  these  applications  strain  relaxation  by  misfit  dislocation 
formation  must  be  avoided.  For  other  applications  it  may  be  possible  to  grow  dislocation  free, 
single  crystal  films  of  difficult-to-grow  semiconductors  through  epitaxial  growth  on  dislocation 
free  substrates  such  as  Si  [5].  Again,  misfit  dislocations  must  be  avoided.  Here  we  review  the 
themoodynamics,  kinetics  and  mechaiusms  of  misfit  dislocation  formation,  paying  particular 
attention  to  the  mechanism  by  which  a  misfit  dislocation  forms  through  the  incremental  movement 
of  a  “threading”  dislocatioa  This  same  mechanism  is  used  to  examine  the  effects  of  an  unstrained 
capping  layer  on  the  velocity  of  a  dislocation  in  a  strained  film.  We  find  that  the  effects  of  an 
unstrained  capping  layer  on  the  dislocation  velocity  are  too  small  to  accoxmt  fw  the  observed 
inhibition  of  misfit  dislocation  formation  in  Si-Oe  films  with  Si  capping  layers.  This  result  leads 
us  to  conclude  that  the  primary  effects  of  such  coping  layers  are  to  inhibit  dislocation  nucleation 
and  multiplication. 

2.  Tlift.  Critical  Thickness  for  Misfit  Dislwation  Formation 

Through  the  pioneering  work  of  van  dcr  Mcrwe  [6]  and  Matthews  [7-9]  it  is  now  well  established 
that  misfit  dislocations  can  fonn  in  heteroepitaxial  films  only  if  the  film  exceeds  a  critical  thickness. 
Here  we  give  a  brief  account  of  the  equilibrium  theory  of  misfit  dislocations;  a  moe  complete 
treatment  of  this  can  be  found  in  a  recent  review  paper  by  Nix  [10].  Consider  a  film  with  an 
equilibrium  lattice  parameter  afiim  deposited  onto  a  substrate  with  a  lattice  parameter  asubstrate 
(^substrate  <  afiim  for  the  Si-Ge  system  discussed  in  this  paper).  Because  of  the  lattice  mismatch, 
the  film  must  be  strained  in  order  to  establish  an  epitaxial  relationship  with  the  substrate.  It  is  now 
well  established  that  if  the  film  is  very  thin,  less  that  the  critical  thickness,  the  stable  equilibrium 
state  is  one  in  which  the  film  is  homogeneously  strained  and  contains  no  misfit  dislocations.  Only 
if  the  critical  thickness  is  exceeded  can  the  free  energy  of  the  film  be  reduced  by  the  presence  of 
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misfit  dislocations.  For  such  thicker  films,  introduction  of  misfit  dislocations  lowers  the  free 
energy  by  reducing  the  homogeneous  clastic  strain  in  the  film.  The  equilibrium  state  can  be 
determined  by  finding  the  misfit  dislocation  spacing  at  which  this 
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Fig.  1 .  Critical  film  thickness  versus  misfit  strain  for  Si-Ge  films  on  SL  The  results  of  Bean  et 
aL  [1 1]  are  compared  with  the  equilibrium  theory. 


reduction  in  fiee  energy  is  just  balanced  by  the  increase  in  free  energy  associated  with  the  energies 
of  the  misfit  dislocations  themsdves.  Then,  the  criticai  thickness,  he,  is  that  film  thickness  for 
which  the  equilibrium  misfit  dislocation  spacing  approaches  infinity.  The  resulting  expression  for 
the  critical  thickness  takes  the  form 


(1) 
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where  M.  and  V  are  isotropic  clastic  constants,  b  is  the  Burgers  vector  of  the  misfit  dislocation,  Mf  is 
the  biaxial  elastic  modulus  of  the  film,  B  is  a  constant  and  e  is  the  misfit  strain,  defined  in  the 
present  discussion  to  be 

^film  ^substrate  ^film  *  substrate  0) 

£  _  - - -  ,  - - - 

substrate  film 


Film  Thickness,  h  (nm) 

Fig.  2.  Misfit  dislocation  spacing  in  Si-Gc  films  grown  to  various  thicknesses  on  Si  [1].  The 
observed  dislocadcxi  spadngs  are  compared  with  the  predictions  of  the  equilibrium 
theory. 

The  equilibrium  model  shows  that  the  critical  thickness  depends  strongly  on  the  misfit  strain  e. 

For  Si-Ge  films  on  Si  substrates,  the  misfit  strain  varies  alnoost  lineariy  with  the  atomic  fraction  of 
Ge  in  the  film.  As  a  result,  the  predicted  critical  thickness  depends  on  the  composition  of  the  film, 
as  shown  in  Fig.  1.  The  figure  also  includes  experimental  measurements  of  the  critical  thickness 
reported  by  Bean  eLaL.[l  1].  We  note  that  for  films  containing  between  15  and  30  atomic  percent 
Ge,  the  actual  film  thicknesses  at  which  misfit  dislocations  are  first  obseivni  are  much  greater  that 
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the  critical  thicknesses  prediaed  by  the  equilibrium  theory.  This  same  discrepancy  is  also  shown 
in  Fig.  2,  where  the  mialit  dislocation  spacings  predicted  by  the  equilibrium  theory  are  compared 
with  observed  misfit  dislocation  spacings  for  Si-Ge  films  grown  on  Si.  Th“sc  data  were  obtained 
for  Si  +  0.2Ge  films  grown  on  Si  by  limited  reaction  processing  [1].  We  note  that  the  observed 
misfit  dislocation  spacings  are  very  much  greater  than  the  predictions  of  the  equililrium  theory, 
especially  below  a  film  thickness  of  about  80  nm.  This  apparent  critical  thickness  is  very  much 
greater  than  the  critical  thickness  of  13.5  nm  prediaed  by  the  equilibrium  theory.  We  believe  the 
discrepancies  shown  in  these  figures  are  due  to  the  sluggishness  of  the  processes  by  which  misfit 
dislocations  form.  The  kinetic  factors  responsible  for  these  effects  are  described  in  the  next 
section. 


3.  Kinetic  Theory  of  Misfit  Dislocation  Formation 
The  Dislocation  Velocity 

The  kinetics  of  forming  misfit  dislocations  can  be  studied  by  considering  the  mechanism  by  which 
a  misfit  dislocation  forms.  Following  the  work  of  Matthews  [7*9],  we  consider  u.at  misfit 
dislocations  are  created  by  the  movement  of  “ihreadbg”  dislocations  that  deposit  misfit  dislocations 
at  the  film/substraie  interface  as  they  move.  We  use  the  quotation  marks  here  and  throughout  this 
paper  to  denote  any  dislocation  that  extends  firom  the  film/substrate  interface  to  the  free  surface  of 
the  film;  the  “threading”  dislocation  need  not  have  origmated  in  the  substrate.  Figure  3  shows  a 
“threading**  dislocation  moving  in  a  film  and  depositing  a  misfit  dislocation  as  it  moves.  The 
kinetics  of  movement  of  the  dislocation  can  be  studied  by  using  the  methods  of  Matthews  [7-9]  and 
Freund  [12-13],  who  showed  that  the  work  done  on  the  moving  dislocation  by  the  stress  in  the 
film  must  be  greater  that  the  work  needed  to  form  the  misfit  dislocation.  For  the  case  of  a  Si-Ge 
film  grown  onto  a  (001)  Si  substrate,  the  wwk  done  by  the  stress  in  the  film  for  a  unit 
displacement  of  the  “threading”  dislocation  is 


where  x  is  the  resolved  shear  stress  and  o  the  biaxial  stress  in  the  film,  b  is  the  Burgers  vector,  h  is 
the  film  thickness  and  the  numerical  factors  reflect  the  orientation  of  the  slip  plane  and  slip  direction 
in  the  (001)  film  For  the  (001)  film  under  consideration,  the  misfit  dislocation  is  a  60  degree 
mixed  dislocation.  The  energy  of  a  unit  length  of  such  a  misfit  dislocation  is 
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where  (if  and  Hs  are  the  shear  mcxiuli  of  the  fihn  and  substrate,  respectively,  p  is  a  constant  taken 
to  be  0.755  and  the  other  terms  have  their  usual  meaning  [10].  Freund  [12-13]  has  shown  that 


"Threading” 

Dislocation  distance 


Fig.  3.  Modtxi  of  a  “threading”  dislocation  which  produces  a  misfit  dislocation  at  the  interface 
between  the  strained  film  and  the  substrate.  Movement  of  the  “threading”  dislocation 
by  a  unit  distance  creates  a  unit  length  of  misfit  dislocation. 

the  critical  film  thickness  can  be  derived  by  setting  Wiayer  =  Wgodeg-  This  is  a  condition  of 
equilibrium,  at  which  point  the  “threading”  dislocation  does  not  move.  For  the  “threading" 
dislocation  to  move  at  a  finite  velocity  and  to  deposit  a  misfit  dislocation  as  it  moves,  Wiayo-  must 
be  greater  that  Wgodeg-  We  define  Wnct  =  Wiaycr  *  W60deg  as  the  net  weak  done  on  the  moving 
dislocation  as  it  moves  a  unit  distance.  We  may  think  of  this  net  work  as  being  done  by  an 
effective  shear  stress  defined  by 

t  ^bh  /4  =  W„,  =W,  -W^.  ,  (5) 

eff  V  2  net  layer  60  deg 

This  effective  stress  drives  the  “tiveading”  dislocation  at  a  finite  velocity,  and  leads  to  the 
dissipation  of  the  net  work  as  heat  Using  cqns.  (3)  and  (4)  the  effective  shear  stress  may  be 
written  in  final  form  as 
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47c{1- v)(Hf +  ^s)  h 


In 


(6) 


(7) 


or 


%ff  ^  ^threshold  •  (8) 

H«eafter,  we  use  the  symbol  A  to  describe  the  coefficient  of  the  logarithmic  term: 

because  it  appears  in  so  many  of  the  equatioos.  We  note  that  the  critical  thickness  can  be  found  by 
setting  Xeff  =  0*  For  tl^c  “threading”  dislocation  to  move  at  a  finite  velocity  the  effective  shear 
stress  must  be  greater  than  zero,  Xeff  >  0. 

The  effective  shear  stress  develqied  here  is  equivalent  to  the  “excess”  stress  disciissed  by  Dodson 
and  Tsao  [14-15].  It  is  that  part  of  the  shear  stress  which  is  not  balanced  by  the  drag  of  the  misfit 
dislocation  and  is  available  to  drive  the  “threading”  dislocation  forward.  We  use  the  effective 
stress  terminology  here  because  it  is  well  established  in  the  literature  on  creep  of  diamond  cubic 
crystals  [16]. 

The  rate  of  motion  of  the  “threading”  dislocation  in  the  Si-Ge  film  is  assumed  to  be  governed  by 
the  thermally  activated  nucleari<»  and  motion  of  Idnks  along  the  dislocation  line.  Because  these  are 
near-atomic  scale  processes,  it  is  expected  that  the  kinetics  of  dislocation  motion  in  thin  film 
structures  will  be  the  same  as  for  bulk  crystals.  Accordingly,  based  on  experimental  work  of 
Alexander  and  Haa^n  [17],  we  assume  that  the  “threading”  dislocation  velocity  can  be  expressed 
as 
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(10) 
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where  U  is  the  activation  energy  for  dislocation  motion,  B  (=  7.33x10^  m/s)  and  to  (=9-8Mpa)  are 
constants  and  teff  is  the  effective  stress  for  dislocatior  motion  given  by  eqn.  (6).  Activation 
energies  for  dislocation  motion  in  Si  and  Ge  were  measured  by  Alexander  and  Haasen  using  etch 
pitting  methods.  Their  results  are  shown  in  Table  L  Also  shown  in  Table  I  are  activation  energies 
for  dislocation  motion  in  Si-Ge  films,  reported  recently  by  Hull  ct  al..  [18]  and  Tuppen  and 
Gibbings  [19].  The  results  of  Hull  et  al..[18]  arc  based  on  in-situ  annealing  of  TEM  samples 
whereas  the  results  of  Tuppen  and  Gibbings  [  19]  arc  based  on  a  scratch/etching  method.  The 
results  of  Hull  et  aL.[18]  suggest  that  the  activation  energy  might  depend  strongly  on  corr^sition 
or  strain.  We  have  also  conducted  in>situ  TEM  annealing  experiments  on  Si  1-xGex  films  with  the 
following  compositions:  x  =  0.16, 0.22, 0.34.  We  find  the  activation  energy  for  motion  to  be  2.2 
±  0.2  eV  for  all  compositions  (and  strains)  studied  Some  of  these  experimental  results  are  shown 
in  Fig.  4. 


Table  I 

Reported  Activation  Energies  for  Dislocation  Motion 
in  Si.  Ge  and  Si-Ge  Epitaxial  Thin  FUms 


Material 

Technioue 

Bulk  Si 

Etch  Pitting 

BulkGe 

Etch  Pitting 

Strained  Sii-xGcx 

Layers  (x=0.3) 

In-SituTEM 

Annealing 

Strained  Sii-xGcx 

Layers  (x=0.05,0.13) 

Scratch/ 

Etching 

Strained  Sii-xGcx 

Layers  (x=0. 16,0.22,0.34) 

In-SituTEM 

Annealing 

Activation 

Reference 

Energy 

2.2  eV 

Alexander  and  Haasen  [17] 

1.6  eV 

Alexander  and  Haasen  [17] 

1.110.2  cV 

Hull  et  al .  [18] 

2.210.2  eV 

Tuppen  and  Gibbings  [19] 

2.210.2  eV 

present  work 

A  few  comments  should  be  made  about  the  discrepancy  that  exists  between  the  results  of  Hull  et 
aL.  [18]  and  the  present  wexit  Two  different  observational  techniques  have  been  used  to  determine 
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the  activation  energies  for  dislocation  motion  in  the  in-situ  TEM  annealing  experiments.  The 
technique  used  by  Hull  et  al..  (18),  and  also  in  preliminary  experiments  in  the  present  work, 
involves  the  obscrvaticwi  of  a  particular  area  of  the  foil  and  the  mcasuretr«nt  of  the  velocities  of  the 
dislocations  that  happen  to  come  into  view.  We  found  an  activation  energy  of  1.3- 1.6  eV  for  a 
film  with  22%  Ce  using  this  technique.  A  diflferent  techiuque  was  used  to  obtain  the  experimental 
data  shown  in  Fig.  4.  In  this  method  a  suitable  dislocation  is  found  when  the  temperature  is  low 
and  its  motion  is  monitored  as  the  .sample  is  heated  to  progressively  higher  temperatures.  Because 
dislocation  velocities  are  so  fast  at  high  temperanires  the  individual  dislocation  measurements  could 
be  carried  out  wily  in  the  range  of  520-67 5*C  Higher  temperature  measurements  were  taken  from 
randomly  captured  events.  This  method  gives  an  activation  energy  of  2.2  eV,  even  for  a  film  with 
30%  Gc. 

The  technique  based  on  the  observation  of  a  particular  area  of  the  foil  may  give  cnoncous 
activation  energies  because  at  lower  temperatures  only  the  fastest  dislocations  are  observed, 
whereas  at  high  temperatures  all  dislocations  in  the  velocity  spectrum  may  be  sampled  This  could 
produce  the  low  apparent  activation  energies  that  have  been  observed 

Because  the  single  dislocation  in-situ  TEM  annealing  experiments  give  results  that  are  in  agreement 
with  the  etch  pitting  results  of  Tuppen  and  Gibbings  [19]  and  Alexander  and  Haasen  [17],  we 
conclude  that  the  activation  energy  for  motion  of  “threading”  dislocations  in  Si-Ge  films  is  about 
2.2  eV. 

EYQluBgaof  ihc..Mi5fiiC>islQcation  Dcnsiiy 

We  may  now  calculate  the  rate  at  which  misfit  dislocations  are  fonned  in  an  epitaxial  Si-Ge  film 
with  a  fixed  density,  N,  of  “threading”  dislocations.  Because  each  “threading”  dislocation  creates 
a  misfit  dislocation  as  it  moves,  it  follows  that  the  mean  spacing,  S,  between  misfit  dislocations 
changes  with  time  according  to 


dt 


(11) 
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Fig.  4.  In-Situ  high  voltage  TEM  measuretoents  of  ‘threading”  dislocation  velocities  in  Si-Ge 
films.  Measurements  are  based  on  the  technique  of  measuring  the  motion  of  an 
individual  dislocation  at  various  temperatures.  The  measured  achvation  energy  is  2.2 
eV. 

where  v  is  the  “threading”  dislocation  velocity  given  by  eqn.  (10).  These  equations  can  be  used  to 
predict  the  evolution  of  the  misfit  dislocation  spacing  during  film  growth.  Results  for  different 
“threading”  dislocation  densities  are  shown  in  Fig.  5  and  conq>ared  with  equilibrium  predictions. 
We  note  that  for  a  low  “threading”dislocation  density  the  predicted  misfit  dislocation  spacings  are 
much  larger  than  the  equilibrium  spacings.  However,  these  predictions  for  a  constant  “threading” 
dislocation  density  are  unlike  the  data  shown  in  Fig.  2.  A  rapid  reduction  in  tl»  misfit  dislocation 
spacing  at  a  film  thickness  of  about  80  nm  is  not  predicted.  We  believe  that  this  rapid  change  in  the 
misfit  dislocation  spacing  is  caused  by  the  nuclcadon  of  new  “threading”  dislocations  and  the 
multiplication  of  existing  ones. 
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Following  methods  developed  in  the  creep  literature  [16]  we  may  assume  that  “threading” 
dislocadons  multiply  according  to  a  law  of  the  form 


mulnpucanon 


(12) 


where  5  is  a  breeding  factor.  When  this  law  is  used  in  conjunction  with  eqn.  (1 1)  the  misfit 
dislocation  spacing  is  predicted  to  change  more  abrupdy  during  film  growth.  Results  for  various 
constant  breeding  factors  are  shown  in  Fig.  6.  We  note  that  a  sharp  reduction  in  the  misfit 
dislocation  spacing  is  prediaed  to  occur  near  a  film  thickness  of  lOOnm. 


Fig.  5.  G^tqputed  misfit  dislocation  spacings  during  the  growth  of  a  Si-Ge  filni  (xi  Si  fes* 

different  constant  '^threading’'  dislocadon  densities.  The  results  are  compared  with  the 
predictions  (Xf  the  equilibrium  theory. 
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Fig.  6.  Confuted  misfit  disIcx»tion  spacings  during  the  growth  of  a  Si-Gc  film  on  Si  for  various 
constant  breeding  factors..  The  results  are  con:q)ared  with  the  predictiOTS  of  the 
equilibrium  theory. 


It  is  now  well  established  diat  misfit  dislocations  can  fonn  during  high  temperature  annealing  of 
films  that  were  initially  dislocation  free.  This  indicates  strcmgly  that  dislocation  niK^leation,  which 
probably  occurs  at  the  free  surface  of  the  film,  also  contributes  significantly  to  the  formation  of 
new  “threading”  dislocations.  Such  dislocation  nucleation,  although  not  well  understood,  is 
expected  to  cause  the  misfit  dislocation  spacing  to  decrease  abruptly  with  increasing  film  thickness. 


Kinetics  of  Strain  Relaxation 

The  processes  of  misfit  dislocation  formation  discussed  here  lead  to  a  relaxation  of  dre 
homogeneous  elastic  strain  in  the  film.  This  relaxation  may  be  expressed  as  a  biaxial  plastic  strain 
rate  of  the  form 
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(13) 


e  ,  .  =|Nbv 

plastic  ^ 


or  a  stress  relaxation  rate  of  the  form 

^  Elastic  ■ 

These  equations  can  be  used  to  describe  the  Idnetics  of  strain  relax^on  at  elevated  temperatures, 
provided  they  are  coupled  with  an  equation  for  the  evolution  of  the  “threading”  dislocation  density. 
Assuming  that  both  nucleation  and  multiplication  of  “tiireading”  dislocations  occurs,  an 
evolutionary  equation  of  the  fcnm 

. 

nucleation 


can  be  used. 


4.  Effects  of  Unstrained  Caopiny  Lavers  on  the  Processes  of  Misfit 

PistotatiQn  Fqcmatiffl 

Unstrained  epitaxial  capping  layers  on  strained  epitaxial  films  are  a  natural  feature  of  some  device 
structures.  Such  layers  are  known  to  inhibit  the  formaticm  of  misfit  dislocations  within  the  strained 
film.  Here  we  craisider  die  physical  processes  by  which  the  formation  of  misfit  dislocations  mi^t 
be  inhibited. 


Effects  of  Unstrained  Capping  Lavers  on  the  Motion  of  •Threading”  Dislocations 

Here  we  cmsider  the  influence  of  unstrained  coping  layers  on  the  motion  of  “threading” 
dislocations  in  strained  films.  In  each  of  the  cases  to  be  discussed,  an  effective  shear  stress  for 
dislocation  naotion  will  be  determined.  These  effective  stresses  will  lead  directly  to  predicted 
velocities  through  the  use  of  eqn.  (10). 
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We  consider  first  an  unstrained  capping  layer  which  is  sufficiently  thin  that  the  “threading” 
dislocation  is  able  to  drag  through  the  cap  as  it  moves.  This  is  depicted  in  Fig.  7.  Because  the 
capping  layer  has  the  same  lattice  parameter  as  the  substrate,  it  is  not  elastically  strained  and  it  does 
not  exert  forces  on  the  pan  of  the  “threading”  dislocation  that  lies  in  the  cs^ing  layer.  As  a 
consequence,  the  pan  of  the  “threading”  dislocation  that  extends  through  tte  capping  layer  must  be 
dragged  along  by  the  pan  of  the  dislocation  that  lies  in  the  strained  film.  As  shown  in  the  figure, 
forces  can  be  transmitted  to  the  dislocation  segment  in  the  unstrained  capping  layer  through  line 
tension  effects.  The  effective  stress  for  dislocation  motion  in  such  a  capped  structure  can  be 
expressed  as 


where  the  notation  “cut”  signifies  that  the  “threading”  dislocation  cuts  through  the  cap,  t  is  the 
thickness  of  die  capping  layer  and  the  other  terms  are  the  same  as  those  defined  for  eqn.  (6).  This 
equation  differs  for  eqn.(6)  in  two  ways.  Hrst,  the  logarithmic  term  includes  the  total  thickness 
the  strained  film  and  the  capping  layer,  hft,  because  the  “threading”  dislocation  moves  through 
both  layers.  Second,  the  term  involving  the  shear  stress  x  is  reduced  by  the  fiaction  h/fb-ft)  to 
reflea  the  faa  that  work  is  done  only  aa  that  part  of  the  “threading”  dislocation  diat  lies  in  the 
strained  film. 

If  the  thickness  of  the  capping  layer  is  greater  than  a  critical  thickness,  then  the  moving  dislocation 
does  not  cut  through  the  cap  and  is  ccHifined  to  move  within  the  strained  film  This  is  illustrated  in 
Fig.  8.  For  this  case,  two  dislocations  are  left  behind  the  moving  dislocation,  tme  at  the 
film/substrate  interface  and  one  at  the  interface  between  the  film  and  tiie  capping  layer.  In  this  case 
the  effective  shear  stress  for  dislocation  motioi  is  defined  in  terms  of  the  work  done  on  the  moving 
dislocation  by  the  stress  in  the  film,  less  the  woik  needed  to  form  the  two  misfit  dislocations: 

■  <■’> 

The  two  identical  logarithmic  terms  are  associated  with  die  energies  of  the  two  misfit  dislocations. 
The  numerical  fturtors  in  die  argument  of  die  logarithmic  terms  are  needed  to  conpute  the  cut  off 
radii  for  the  two  dislocations  in  the  dipole. 
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Fig.  7.  Motion  of  a  “threading”  dislocation  that  cuts  through  a  thin  capping  layer.  The  capping 

layer  is,  like  the  substrate,  unstrained.  Forces  on  the  segment  that  extends  through  the 
c^ing  layer  are  provided  by  line  tension  effects.  The  velocity  of  the  dislocation  is 
v(cut). 


Fig.  8.  Motion  of  a  threading  that  is  confined  to  remain  within  the  strained  film.  The  capping 
layer  is  too  thick  for  the  dislocation  to  cut  through.  The  velocity  of  the  dislocation  is 
v(dipole). 
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Finally,  we  consider  the  rate  of  motion  of  a  “threading”  dislocation  that  extends  fix)m  the 
film/capping  layer  interface  to  the  free  surface  of  the  capping  layer,  as  shown  in  Fig.  9.  The 
motion  of  the  dislocation  segment  is  not  driven  by  the  stress  in  the  film.  Rather  it  is  driven  by  the 
differences  in  energy  of  the  misfit  dipole  and  the  single  misfit  dislocation  left  at  the  film/substrate 
interface  when  the  segment  in  question  moves  by.  For  this  kind  of  movement  and  for  the  case  in 
which  the  capping  layer  is  thin  compared  to  the  film,  the  effective  shear  stress  can  be  expressed  as 

.  (.8) 

where  the  term  “catch-up”  refers  to  die  fact  that  the  segment  in  question  is  attempting  to  catch-up 
v^ith  the  segment  that  deposits  the  dipole  in  its  modem.  The  first  two  terms  in  the  Ixackets  relate  to 
the  energies  of  the  two  dislocations  in  the  dipole  while  the  last  term  relates  to  die  energy  of  the 
single  dislocation  left  behind  by  the  moving  dislocation  segment 


catch-up  with  the  leading  edge  of  the  dipole  not  shown  in  the  figure.  The  velocity  of 
the  dislocation  is  v(catch-up). 


Equations  (16),  (17)  and  (18)  have  been  used  in  conjunction  with  eqn.  (10)  to  examine  the 
influence  of  capping  layers  on  the  velocities  of  the  various  dislocation  segments  under 
consideration.  One  set  of  results  is  shown  in  Fig.  10  where  the  different  velocities  are  shown  as  a 
function  of  capping  layer  thickness  for  a  50  nm  thick  Si-Ge  film.  The  velocities  are  normalized  in 


15 


such  a  way  that  the  velocity  for  an  uncapped  film  is  0.64,  as  shown  by  the  arrow.  Consider  fust 
the  velocity  v(cut)  associated  with  the  motion  of  a  “threading”  dislocation  that  cuts  through  the 
capping  layer.  This  velocity  falls  continuously  with  increasing  capping  layer  thickness.  When  the 
capping  layer  is  3  or  4  times  thicker  than  the  film,  v(cut)  is  about  an  order  of  magnitude  smaller 
than  the  velocity  far  an  uncapped  film.  However,  for  capping  layers  of  this  thickness,  the 
dislocation  is  not  able  to  move  by  cutting  through  the  capping  layer.  Rather  the  segment  that  lies  in 
the  strained  film  will  move  more  quickly  than  the  segment  that  lies  in  the  capping  layer,  with  the 
result  that  a  dipole  configuration  is  eventually  formed.  Figure  10  also  shows  calculated  velocities 
for  the  dipole  configuration.  The  velocity  v(dipoIe)  decreases  with  increasing  capping  layer 
thickness  and  is  smaller  than  v(cut)  up  to  a  critical  capping  layer  thickness  of  about  35  nm. 

Beyond  that  capping  layer  thickness,  vfdipole)  is  constant  because  the  energy  of  the  dipole  does 
not  depend  on  the  thickness  of  the  capping  layer.  The  variation  of  v(dipoIe)  with  capping  layer 
thickness  at  small  capping  layer  thicknesses  is  due  to  the  fact  that  the  dipole  energy  depends  on 
capping  layer  thickness  in  that  regime.  A  comparison  of  v(cut)  and  v(dipole)  indicates  that  cutting 
is  the  dominant  process  at  small  capping  layer  thicknesses,  as  expected,  and  that  dipole  formation 
dominates  at  large  capping  layer  thicknesses.  The  overall  result  is  that  the  capp'-'g  layer  decreases 
the  dislocation  velocity  by,  at  most,  a  factor  of  2. 

It  is  of  interest  compare  the  velocity  v(calch-up)  with  v(cut)  and  vfdipole).  Rgure  10  shows  that 
below  the  critical  capping  layer  thickness  of  35  nm,  v(catch-qp)  is  greater  that  both  v(cut)  and 
v(dipole),  whereas  at  thicknesses  greater  that  the  critical  thickness,  v(catch-up)  <  v(cut)  < 
v(dipole).  This  result  is  understandable  and  is  to  be  expected.  If  a  dipole  configuration  existed  in 
a  sample  with  a  very  thin  ct^iping  layer,  then  the  segment  in  the  capping  layer  would  quickly  catch¬ 
up  with  the  leading  edge  of  the  dipole  and  convert  die  dislocation  into  a  cutting  type  “threading” 
dislocation. 

A  final  note  should  be  made  about  the  present  analysis.  The  calculations  we  have  made  ^ply  to 
strained  films  that  have  not  yet  relaxed.  In  this  limit  the  capping  layer  is  under  no  biaxial  sness. 
During  the  course  of  strain  relaxation  in  films  with  thin  capping  layers,  the  moving  dislocations 
plastically  deform  both  the  film  and  the  capping  layer.  This  relieves  the  strain  in  the  strained  layer 
but  it  causes  a  strain  of  opposite  sense  to  develop  in  the  capping  layer.  This  effect,  which  was  not 
taken  into  account  in  the  above  analysis,  would  reduce  the  critical  capping  layer  thickness  for 
forming  the  dipole  configuration  and  would  slow  the  speed  of  cutting  type  dislocations. 
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Capping  Layer  Thickness  (nm) 

Fig.  10  A  comparison  of  various  dislocation  velocities  as  a  function  of  capping  layer  tiiickncss. 

The  critical  capping  layer  thickness  above  which  dislocations  move  by  the  dipole 
process  is  about  35  nm. 


The  analysis  described  in  the  previous  section  indicates  that  “threading”  dislocaticm  velocities  are 
reduced  by  only  as  much  as  a  factor  of  2.  As  we  shall  see,  this  small  effect  cannot  account  for  the 
large  inhibiting  effects  of  unstrained  capping  layers  on  the  formation  of  misfit  dislocations. 


A  full  accounting  of  the  effects  of  unstrained  capping  layers  must  include  the  effects  on  nuclcarion 
and  multiplication  of  “threading”  dislocations.  Unstrained  Si  capping  layers  of  various  thicknesses 
were  grown  on  dislocation  free  Si-Ge  films  oy  the  LRP  technique.  These  samples  were  then 
annealed  at  850'C  for  4  minutes,  either  in  a  rapid  thermal  oven  outside  of  the  growth  chamber  (ex- 
situ)  or  in  the  LRP  growth  chamber  prior  to  removal  of  the  substrate  (m-situ).  The  misfit 
dislocation  spacings  in  both  capped  and  uncapped  films  were  measured  after  annealing.  The 
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Fig.  1 1.  Effect  of  capping  layer  thickness  on  the  misfit  dislocation  spacing  in  Si-Ge  films  after 
annealing  at  850*C  for  4  minutes.  The  capping  layers  greatly  inhibit  misfit  dislocation 
formation. 

results  are  shown  in  Fig.  1 1.  As  shown  in  the  figure,  the  annealing  process  allows  a  high  density 
of  misfit  dislocations  to  form.  The  misfit  dislocation  spacing  for  the  uncapped  film  is  about  O.S 
pm.  As  can  be  seen  in  the  figure,  the  misfit  dislocation  spacing  after  annealing  increases  markedly 
with  increasing  capping  layer  thickness.  For  a  capping  layer  thickness  of  50  nm,  the  misfit 
dislocation  spacing  is  more  than  100  times  greater  than  the  misfit  dislocation  spacing  fia*  uncapped 
films.  By  comparison,  a  SO  nm  capping  layer  deaeases  the  “threading”  dislocation  velocity  by 
only  about  a  factor  of  2.  This  indicates  strongly  that  something  other  that  the  “threading” 
dislocation  velocity  must  be  affected  by  the  capping  layer.  We  believe  the  primary  effect  of  the 
capping  layer  is  to  inhilnt  misfit  dislocation  nucleation  and  multiplication. 
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B.  In  Situ  Study  of  Isothermal  Strain  Relaxation  in  Si-Ge 
Heteroepitaxial  Films  using  Substrate  Curvature 
Measurements 

Veronique  T.  GUlaid,  David  B.  Noble*  and  William  D.  Nix* 

1.  Introduction 

Lattice  mismatch  between  a  heteroepitaxial  Si-Ge  film  and  a  Si  substrate  causes  elastic  misfit  strain 
to  exist  in  the  film.  When  the  film  thickness  is  higher  than  a  critical  value  [1],  the  biaxial  stress 
associated  with  diis  elastic  strain  can  serve  as  driving  feroe  for  the  movement  of  threading 
dislocations.  As  they  move  through  the  film,  these  threading  dislocations  lay  down  misfit 
dislocations  at  the  film/substxate  interface  and  this  causes  the  elastic  strain  to  relax.  Several  device 
applications  of  Si-Gc  heteroepitaxial  films  are  based  oa  the  significant  reducticar  of  die  indirea 
bandgap  caused  by  the  misfit  strain  [2].  For  these  applications,  relaxation  of  die  misfit  strain  is  not 
desirable.  There  exist  other  applications  for  which  there  is  a  need  for  a  high  density  of  misfit 
dislocations  in  conjunction  with  a  low  density  of  threading  dislocations  [3].  In  both  cases, 
understanding  of  die  kinetics  of  relaxation,  as  well  as  the  mechanisms  involved  is  important  to  die 
development  of  reliable  devices. 

The  biaxial  stress  in  the  film  causes  elastic  beading  of  the  substrate.  The  resulting  substrate 
curvature  is  related  to  the  elastic  strain  in  the  film  through  the  following  equation. 

£=.^J!iLak 

Mf  6hf 

where  M  represents  biaxial  modulus,  h  represents  thickness  and  AK  is  the  difference  in  curvature 
between  bare  substrate  and  substrate  with  fihn.  The  subscripts  “s”and  “f  ”  refer  to  substrate  and 
film,  respectively.  During  isothermal  annealing  experiments,  the  change  in  substrate  curvature 
with  time  is  measured  with  a  laser  scanning  system  [4].  The  relaxation  of  elastic  strain  and  the 
corresponding  evoluticn  of  plastic  strain  are  found  using  the  above  formula. 

The  kinetics  of  fonnation  of  misfit  dislocations  can  be  studied  by  considering  that  they  are 
deposited  at  the  film/substrate  interface  by  threading  dislocations  moving  through  the  film.  Since 
each  threading  dislocation  creates  a  misfit  dislocation  as  it  moves,  the  mean  spacing,  S,  between 
misfit  dislocations  varies  with  time  according  to  the  following  formula  [1]. 
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(2) 


In  the  above  equation,  v  is  the  threading  dislocation  velocity.  N.  the  mobile  threading  dislocation 
density,  and  the  factor  of  2  is  related  to  the  60®  misfit  dislocations.  This  equ^on  shows  that  the 
kinetics  of  misfit  dislocation  formation  depends  on  the  threading  dislocation  velocity  and  the 
density  of  mobile  threading  dislocations,  both  of  which  vary  during  the  course  of  relaxation.  It  is 
necessary  to  separate  these  two  parameters  and  their  evolution  with  time  in  order  to  better 
understand  the  kinetics  of  relaxation. 


The  kinetics  of  movement  of  the  threading  dislocations  can  be  studied  using  an  argument  which 
considers  that  the  work  done  on  the  moving  dislocation  by  the  stress  in  the  film  has  to  be  greater 
than  the  work  necessary  for  the  formation  of  the  misfit  dislocation.  This  method  leads  to  a  value 
for  the  effective  shear  stress,  which  is  the  stress  still  available  to  drive  the  threading  dislocation 
forward  after  a  misfit  dislocation  has  been  deposited  in  its  wake  [1]: 
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where  b  is  the  Burgers  vector  of  the  misfit  dislocation,  v  is  Poisson’s  ratio,  m  and  pf  are  the  shear 
moduli  of  the  substrate  and  film,  respectively,  and  P  is  a  constant  equal  to  0.755  [1]. 
o  is  the  biaxial  stress  in  the  film,  x  is  the  resolved  applied  shear  stress  and  the  expression  appearing 
after  the  minus  sign  represents  the  pan  of  the  shear  stress  which  is  required  for  the  formation  of  the 
misfit  dislocation.  Since  the  vaiiatitm  of  the  elastic  strain  with  time  during  the  annealing 
experiment  is  known,  the  corresponding  variation  of  the  resolved  shear  stress  t  is  also  known  and 
this  allows  for  a  determination  of  the  effective  shear  stress  as  a  function  of  time  using  equation 
(3). 

If  the  rate  of  motion  of  threading  dislocations  is  assumed  to  be  governed  by  the  thermally  activated 
nucleation  and  motion  of  kinks  along  the  dislocation  line,  the  empirical  equation  developed  by 
Alexander  and  Hassen  [5]  may  be  used  to  relate  the  threading  dislocation  velocity  to  the  effective 
stress; 


v  = 


(4) 
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U  is  the  activation  energy  for  dislocation  motion,  n,  B  and  xq  arc  constants,  tejff  is  the  effective 
stress  given  by  equation  (3)  and  the  other  terms  have  their  usual  meaning.  In  the  present  study,  the 
activation  energy  and  the  constants  are  determined  in  films  similar  to  those  studied  in  the  annealing 
experiments  using  a  least  square  method  to  fit  velocities  measui«i  by  TEM  and  effective  shear 
stresses.  Since  the  nucleation  and  motion  of  kinks  is  a  near-atomic  scale  process,  the  use  of 
equation  (4),  which  was  developed  for  bulk  material,  is  expected  to  be  applicable  for  thin  films  [1]. 
Furthermore,  as  will  be  shown  later  in  this  paper,  this  equation  gives  a  good  mathematical 
description  of  our  experimental  data. 

As  mentioned  earlier,  the  formation  of  misfit  dislocations  causes  the  homogei^us  elastic  strain  in 
the  film  to  relax.  The  Oiowan  equation  may  be  used  to  express  the  biaxial  plastic  strain  rate  in  the 
following  manner 


1  4  £ 

e  =  lNvb  =>  N  =  ~ 

4  V  b  (5) 

Since  the  variation  of  the  effective  shear  stress  with  time  during  the  annealing  experiment  is  known 
(equation  (3)),  the  conespcxiding  variation  of  the  threading  dislocation  velocity,  v,  with  time  is 
determined  using  equation  (4).  The  evolution  of  the  strain  rate,  e,  is  also  determined  fiom  the 
evoluti  on  of  strain,  e,  as  a  fimetion  of  time.  Using  equation  (5)  then  allows  us  to  obtain  the 
evolution  of  N,  the  mobile  dislocation  density,  during  the  course  of  annealing. 

There  are  several  physical  mechanisms  which  can  can.se  the  mobile  dislocation  density  to  change 
during  the  relaxation  process.  N  can  increase  if  threading  dislocations  are  nucleated  and/or  if 
dislocation  multiplication  is  occurring;  these  phenomena  are  predominant  in  the  eariy  stage  of 
relaxation  when  the  available  driving  force  ts  still  high.  The  mobile  dislocation  density  can 
decrease  if  threading  dislocations  are  lost  at  the  edge  of  the  film  and/or  if  dislocaticm  interactions 
cause  immobilization.  The  shape  of  the  curve  N(t)  is  determined  by  the  mechanism(s)  causing  the 
mobile  dislocation  density  to  change.  For  example,  at  the  start  of  the  relaxation,  the  general  shape 
of  N(t)  depends  on  whetiier  nucleation  or  multiplication  of  dislocation  is  the  dominant  process,  as 
shown  in  Fig.  1. 
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Nucleation 


Multiplication 


Hgure  1.  General  shape  of  mobile  dislocation  density  versus  tune  with  nucleation  or 
multiplication  being  the  predominant  mechanism  of  dislocation  creation. 

If  nucleation  is  the  dominant  process  for  threading  dislocation  formation,  the  mobile  dislocation 
density  should  start  to  increase  immediately  and  continue  to  increase  at  a  fairly  constant  rate.  On 
the  other  hand,  if  multiplication,  through  a  mechanism  such  as  Hagen-SminJc  [6],  doninates,  it 
will  take  a  certain  amount  of  time  before  a  source  is  created  and  N(t)  will  remain  flat  until 
multiplication  can  start,  then,  N  is  expeaed  to  ircrease  exponentially  with  time  [1], 

=N«v 


Here  5  is  a  breeding  factor  as  in  the  creep  literature.  The  ability  to  determine  the  evolution  of  the 
mobile  dislocation  density  with  time  during  armealing  experiments  should  therefore  help  to  develop 
a  better  understanding  of  the  mechanisms  of  relaxation. 


2.  Sample  and  Experiment  Characteristics 


The  Si-Ge  films  were  grown  in  Prof.  Gibbons’  laboratory  at  Stanford  on  (001),  heavily  n  doped, 
3”  Si  wafers.  An  undoped  Si  layer  was  used  as  a  buffer  between  the  Si-Ge  film  and  the  doped 
substrate.  The  growth  was  performed  at  625  '’C  in  a  lamp-heated  CTVD  reactor  in  the  reaction- 
limited  regime.  The  6  samples  used  in  this  study  have  a  Ge  fraction  of  34.5%  (±  0.5%)  and  film 
thicknesses  of  100, 400, 500, 600, 700,  and  800  A.  Both  Ge  content  and  film  thickness  were 
measured  after  growth  by  RBS.  The  100  A  film  thickness  was  verified  by  cross-section  TEM. 
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After  growth,  a  1  by  1  inch  piece  was  cut  from  the  center  of  the  wafer  to  ensure  film  thickness 
uniformity  over  the  sample.  The  substrate  was  then  thinned  to  about  200  pm  in  order  to  increase 
the  ctuvanire  and  thereby  improve  the  strain  measurement  resolution.  TEM  observations  of  the 
samples  showed  that  all  films,  except  the  100  A  thick  film,  contain  a  certain  number  of  misfit 
dislocations  after  film  growth.  The  100  A  film  is  very  near  the  critical  thickness  for  relaxation. 

The  elastic  strain  in  the  film  after  growth  was  determined  using  a  combination  of  x-ray  diffraction 
and  wafer  curvature  measurements. 

An  annealing  temperature  of  700  *C  was  chosen  and  the  annealing  procedure  was  as  follows:  a 
sample  was  inserted  in  the  furnace  previously  stabilized  at  700  ®C  in  an  argon  atmosphere  and 
reached  the  target  temperature  within  2  minutes.  Starting  immediately  after  insertion,  'ubstrate 
curvature  was  measured  in  situ  every  2  minutes  for  several  hours.  Readers  should  refer  to  von 
Preissig,  [4]  for  a  detailed  description  of  the  furnace  and  the  curvature  measurement  system 

3.  Ejtpgriniaital  Results 

A  calculation  of  the  misfit  strain,  em.  in  a  Si-Ge  film  containing  34.5  %  Ge,  gives  a  compressive 
biaxial  strain  of  -1.30  %  at  room  temperature  and  -1.37  %  at  700  ®C  These  results  are  based  on  a 
film  lattice  parameter  which  takes  into  account  the  negative  deviation  from  Vegard’s  law  [7].  The 
misfit  strain  at  700  X  is  determined  using  experimental  data  for  the  variation  of  the  linear 
expansiwi  coefficient  of  Si  and  Ge  with  temperature  [8-9]  and  a  simple  rule  of  mixtures  for  the 
linear  expansion  coefficient  of  Si-Ge.  Using  a  direct  measurement  of  the  Si-Ge  lattice  parameter  at 
700  °C  for  a  Ge  fraction  of  34.7%  [7]  leads  to  a  value  of  - 1.36  %. 

Experimental  values  for  the  misfit  strain  and  the  elastic  strain  remaining  in  the  films  after  growth 
are  presented  in  Table  I. 

Measured  values  of  strain  are  not  available  for  the  100  A  film  One  can  see  in  Table  I  that  the 
experimental  value  obtained  for  the  misfit  strain  at  room  temperature  ( Em  R.T. )  is  very  close  to 
the  calculated  value  given  above  (-1.30  % ).  However,  the  misfit  strain  measured  at  7(X)  °C,  ( Em 
700°C )  is  significantly  smaller  than  that  calculated  on  the  basis  of  lattice  parameter  measurements 
published  in  the  literature  [7],  ( - 1.36  %  ).  This  effect  does  not  seem  to  be  caused  by  the  curvature 
measurement  system,  as  the  change  in  curvature  due  to  temperature  gradients  through  the  substrate 
thickness  has  been  previously  measured  [4]  and  leads  to  strain  variations  which  are  an  order  of 
magnitude  smaller  than  those  observed  here.  Furthermore  this  effect  is  not  observed  on  a  bare  Si 
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substrate.  Therefore,  it  seems  that  as  the  samples  are  heated  to  700  ®C,  a  small  part  of  the  clastic 
strain  is  immediately  relaxed  even  in  the  case  of  a  film  which  exhibits  no  further  relaxation  for 
about  an  hour,  ( 400  A  film )  as  shown  later.  We  are  unable  to  explain  this  phenomenon  at 
present  Sandies  similar  to  those  described  in  this  study  but  grown  by  MBE  show  the  same  type 
of  effect 


Table  I 

Experimental  values  of  strain  [%]. 


hf  [A] 

ecL(l)  R-T. 

ed.(l)  700®C 

£pl 

Em  R-T. 

Em  700®C 

400 

-1.26 

-1  21 

-0.00 

-1.26 

-1.21 

500 

-1.25 

-1.19 

-0.00 

-1.25 

-1.19 

600 

-1.25 

-1.21 

-0.01 

-1.26 

-1.22 

700 

-1.28 

-1.22 

-0.04 

-1.32 

-1.26 

800 

-1.12 

-1.07 

-0.13 

-1.25 

-1.20 

Eel.(l)  R.T.  is  the  remaining  elastic  strain  after  growth  measured  with  a  combination  of  x- 
ray  diffiaction  method  and  wafer  curvature  at  room  tecqperatuie. 

EeL(l)  700*’C  is  the  elastic  strain  measured  with  the  wafer  curvature  technique  immediately 
after  sample  insertion  in  the  furnace.  Sample  tonperatuie  is  already  within  S 
of  the  target  temperature  at  the  time  of  measuremem. 

£pl  is  the  plastic  strain  CDtresponding  to  the  amount  of  relaxaticMi  which  todc  place  during 
growth  determined  using  TEM  measurements  of  average  misfit  dislocation  spacing,  S, 
(Epi  =  b/2S). 

(Em  R.T.) » (eei.(i)  R.T.)  +  (Epi )  and  (em  7(X)®C )  =  (eci.(i)  700®C )  +  (Ej^  ) 


The  results  of  the  annealing  for  the  6  samples  are  presented  in  Fig.  2.  This  figure  shows  that,  as 
expected,  the  100  A  film,  which  is  very  near  the  critical  thickness,  does  not  substantially  relax. 

The  400  A  film  shows  a  l<mg  plateau  with  no  significant  relaxation,  then,  the  elastic  strain  starts  to 
decrease  very  rapidly  until  relaxation  decelerates  and  a  new  phueau  is  reached.  The  500  A  film 
behaves  in  a  very  similar  marmer  but  the  relaxation  occurs  more  quickly.  This  is  expected  since  the 


25 


strain 


driving  force  is  higher  in  this  thicker  film.  In  the  other  films,  no  plateau  is  observed  as  strain 
relaxation  had  already  started  during  growth. 


Figure  2.  Residual  elastic  strain  versus  time  during  annealing  at  700  °C. 


In  the  case  of  the  800  A  film,  most  of  the  relaxation  already  occurred  before  annealing  and  the 
curve  recorded  only  shows  the  decelerating  part  of  the  strain  relaxation.  The  variation  of  the  shape 
of  the  relaxation  curves  with  film  thickness  is  consistent  with  the  fact  that  the  driving  force  for 
strain  relaxation  increases  as  the  film  becomes  thicker. 

Figure  3  shows  experimental  data  of  dislocation  velocities  versus  effective  stress  for  bulk  material 
[10]  and  for  thin  films  (TEM)>  Both  sets  of  data  lie  on  the  same  line,  which  represents  the  best 
least  square  fit  to  equation  (4).  This  shows  that  the  equation  developed  by  Alexander  and  Hassen 
gives  a  good  representation  of  our  velocity  data.  For  the  following  part  of  this  study,  we 
determine  the  constants  in  equation  (4)  by  fitting  Noble’s  data  cwJy.  since  they  are  measured  in 
films  very  similaT  to  those  studied  here,  which  contain  a  certain  amount  of  oxygen.  We  get  the 
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following  values:  the  activation  energy  U  is  2.1  eV,  the  exponent  n  is  0.86,  and  the  other  constants 
are  B  =  19300  m/s  and  to  =  9  8  MPa. 


Rgure  3.  Variation  of  dislocation  velocity  with  effective  stress. 


Using  the  measurement  of  strain  versus  time,  we  detenxune  the  variation  of  strain  rate  during  the 
course  of  relaxation.  This  data  is  presented  in  Hg.  4  for  the  400  A  thick 
film.  Also  presented  in  this  figure  is  the  variation  of  dislocation  velocity  with  annealing  time 
determined  using  equation  (4)  with  the  values  given  above  for  the  constants.  Then,  the  mobile 
dislocation  density,  N,  is  determined  with  equadem  (5)  and  leads  to  the  curves  displayed  in  Fig.  5. 
In  this  figure,  we  see  that  for  the  thickest  film,  where  most  relaxadon  took  place  during  growth, 
the  nx>bile  dislocation  density  remained  rather  constant  during  annealing.  The  thinner  films  show  a 
sharp  increase  in  N  until  a  peak  is  reached,  after  which  the  mobile  dislocation  denary  decreases 
very  rapidly  to  reach  a  plateau.  The  shape  of  the  curve  at  the  stage  where  N  increases  suggests  a 
multiplication  mechanism.  The  decrease  in  N  should  be  due  to  a  cennbination  of  dislocation 
annihilation  at  the  film’s  edges  and  immobilization  due  to  interaction  with  other  dislocations. 

Using  samples  with  different  areas  should  allow  us  to  determine  which  of  these  mechanisms  is 
dominanL 
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Figure  4.  Strain  rate  and  dislocation  Figures.  Mobile  dislocation  density  versus 

velocity  versus  time  for  the  400  A  film.  time  for  the  400, 500  and  800  A  films. 


4.  Conclusions 

The  combination  of  curvature  measurements  taken  during  annealing  experiments  and  previous 
TEM  dislocation  velocity  measurements  allows  us  to  determine  the  evolution  of  mobile  threading 
dislocation  density  during  the  course  of  strain  relaxation  in  Si-Ge  films.  This  should  prove  to  be  a 
useful  tool  fCR-  a  better  understanding  of  the  kinetics  of  strain  relaxation  in  heteroepitaxial  fixns  and 
the  corresponding  mechanisms  responsible  for  diat  relaxation. 
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